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Abstract 



In this paper, we give a precise description of the rescahng behaviors 
of global strong polynomial solutions to the reformulation of zero sur- 
face tension Hele-Shaw problem driven by injection, the Polubarinova- 
Galin equation, in terms of Richardson complex moments. From past 
results, we know that this set of solutions is large. This method can 
also be applied to zero surface tension Stokes flow driven by injection 
and a rescaling behavior is given in terms of many conserved quantities 
as well. 

Keywords: Hele-Shaw flow, rescaling behavior, polynomial solutions, Richard- 
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1 Introduction 

This present paper is mainly devoted to two problems: two-dimensional zero 
surface tension Hele-Shaw flows driven by injection at the origin and two- 
dimensional zero surface tension Stokes flow driven by injection at the origin. 
In both cases, ingenious complex variable methods have been developed 
which use conformal mapping. Let be the moving domain and f{(,,t) : 
D n{t) is univalent and analytic in D, f{0,t) = 0, / (0, t) > . We set 



For the former problem, the differential equation with injection strength 27r 
at the origin is, 
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f{C,t) = f\tt), D = {x\\x\<l}. 




Equation (1.1) is called the Polubarinova-Galin equation since Galin and 
Polubarinova-Kochina first derived it and investigated the Riemann map- 
ping method along these hnes. A solution to equation (1.1) is said to be 
a strong solution for t G [0,6) if /(^,t) is univalent and analytic in D, 
f(0,t) = 0, /'(C,0) > and continously difFerentiable in t, t e [0,6). For 
the later problem, the differential equation with injection strength 27r at the 
origin is, 

+ 2A'(e,i) = 0, C&dD (1.2) 

where A{^,t) = x(/(C,*)) and x{z) ^ - \ogz as z ^ 0. 
Define 

0{E) = |/ I / is univalent and analytic in E, /(O) = 0, /'(O) > o| . 

In the case of ZST Hele-Shaw flows driven by injection, it is shown in 
Gustafsson [2] that starting with a degree n polynomial mapping /n(^,0) G 
0{D), the strong solution to (1.1) /„(^, t) is also a polynomial of the same de- 
gree. In Huntingford [4] , it is shown that not any given polynomial function 
/(C, 0) G 0(D) can produce a global strong polynomial solution. However, it 
is proven in Gustafsson, Prokhorov and Vasil'ev [3] that starting with a star- 
like mapping /(^, 0) G 0(D), the strong solution to (1.1) f($,,t) is global. In 
the case of ZST Stokes flow driven by injection, in Cummings, Howison and 
King [1], the authors obtain many polynomial solutions to (1.2). The global 
dynamics is not as clear as that of Hele-Shaw flows but some of polynomial 
solutions to (1.2) are global. 

In this paper, we focus on large-time rescaling behaviors of global poly- 
nomial solutions to ZST Hele-Shw flows driven by injection and ZST Stokes 
flow driven by injection in terms of many invariant quantities. These in- 
variant quantities are defined in each problem respectively as follows. Let 
a family of domains {^(t)}t>o be a solution to ZST Hele-Shaw flows driven 
by injection with strength 27r or ZST Stokes flow driven by injection with 
strength 27r and /(^, t) G 0(D) be the corresponding Riemann mapping. In 
the former problem, Richardson complex moments are defined as 

Mk(t) = -[ z^dxdy = ^j f\^,t)f'(i,t)J(-t)d^, z = x + iy. 
The rate of change of moments is given by 

j^Mk(t) = 26ko, 
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d_ 
dt 



where Sko = 1 when A; = and Sko = if A; 7^ 0. In particular, the moments 

Mk, k > 1 are conserved. The zero moment Mo(t) equals ^ | fl{t) \ where 
I i}{t) I is the area of the domain Q(t). In the later problem, it is shown in 
Cummings, Howison and King [1] that the quantities Cfc (t) , /c > where 



1 



^.1 



l«l=i 



satisfy 



2(5, 



In particular, the quantities, Ck,k> 1 are conserved. The quantity Co(t) = 
Co(0) + 2i equals ^ \ n{t) \. 

This current paper is organized as follows. In section 2, we assume 
fn{^,t) = Y!d=i C'ii't)^ 8- global strong degree n > 2 polynomial solution to 
(1.1) and obtain how each coefficient decays or grows in terms of Richardson 
complex moments. By this result, wc obtain many rescaling behaviors for 
fniCit). Let no = min{A; > l|Mfc / 0}, then it is shown that 



lim 

t— »oo L 



and 



lim 

t— >CXD 



fn{V2ic, t) {2t + ^) e - j^MTie 



-M|(0) 



k=2 
k=2 



(V2t) 



Thus, the rescaled limit depends on the initial data only through the con- 
served moments. A geometric characterization of a rescaling behavior for the 
corresponding moving domain = /„(£), i) is given as well. By rescaling 
Q,{t) to be Q'{t) which has area tt only, we obtain that the radius and curva- 
ture of f2'(t) decay to 1 algebraically as ^-("■o/2+i)_ jj^ section 3, it is shown 
that for a global degree n > 2 polynomial solution fn{C,t) = Y17=i to 
(1-2), 



lim 

t— *oo 



/„({,«) -\/2t + Oo(0)4 



E 

k=2 



A geometric characterization of a rescaling behavior for the corresponding 
moving domain = fn{D,t) is given as well. By rescaling to be 
Q,'{t) which has area vr only, we obtain that the radius and curvature of 
Q'{t) decay to 1 algebraically as t~^. 
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2 Rescaling behavior in Hele-Shaw flows 



If we assume that fn{(.,t) = Yl7=i'^i^^ ^ strong degree n polynomial 
solution to the P-G equation (1.1), and that i}{t) = fn{D,t), then by the 
result in Richardson [6], the moments {Mfc(t)}o<fe<n-i can be represented 
by these coefhcients {ak}k>i as 

h,--- ,ik+i 

The zero moment Mo{t) = X^ILi ^ I P= 2i + Mo(0) = al{t) + g{t) where 
g{t) = Y17=2^\^i(^)\^- ^^^^ fn{^,t) is global, by Kuznetsova [5], 

af (t) — 2t is nondecreasing and g{t) is not increasing. Therefore, the coeffi- 
cient ai{t) satisfies 

lim^ = l 



and ak{t),k > 2 is bounded by \/giO) in this case. 

Now given a rcal-cocfficicnt global strong solution ai(t)^ + a2(i)^^ + 
0'3{t)C'^ to (1.1), then by (2.1), a2it) and a-i{t) are expressed in terms of 
ai{t) and the conserved moments {Mi;}i<k<2 as follows: 

a2(t) = Mi- 



-3M2 + af{t) 
and 

The moments in this case are real-valued as well since {afc(t)}i<fc<3 are 
real-valued. We can see that the coefficients {ak{t)}2<k<3 in this example 
satisfy 



and 



Moreover, 



lim a2{t){ai{t)y = Mi = Mi 

t—^oo 



lim a3{t){ai{t)f = M2 = M2. 

t—*oo 



a2it)al{t) - Ml = a2{t)alit) - Mi = - ^^^ ^ ' = O 



3M1M2 
3M2 + afit) 



a3{t)al{t) - M2 = a3{t)al{t) - M2 = 0. 
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We give a general statement of this kind of decay for coefficients of 
global strong polynomial solutions to (1.1) in subsection 2.1, and use the 
results to obtain rescaling behaviors of global strong polynomial solutions 
in subsection 2.2. If the lower Richardson moments vanish, the decay rate 
of coefficients is much faster and the global strong polynomial solutions 
to (1.1) have better rescaling behaviors. In subsection 2.3, a geometric 
characterization of a rescaling behavior for the moving domain is given. 

2.1 Decay of coefficients in terms of moments 

Lemma 2.1. If fn{£,, t) = Y17=i ^ global strong polynomial solution 

of degree n> 2 to the P-G equation (1.1), then for k >2 

lim aiUk = Mfc_i. (2.2) 

t— >oo 

Furthermore, for k >2, 

aW-Wr,=0(^-^). (2.3) 

Proof. We want to prove (2.2) by induction. 

Stepl: Since a„ = Mn-ia^"", it is clear that (2.2) holds for k = n. 

Step2: Assume that the (2.2) holds for n — sq < k < n where sq + I <n — 2. 

Claim: The result (2.2) holds for A; = n - {sq + 1). 

Proof, (of claim) For k = n — (sq + 1), 

Mk-i{t) = ^ iitti^Oi^ ■ ■ ■ ai^ai^+...+i^ 

ii,--- ,ik 

=a\ak+ iiOi^ai^ ■ ■ ■ Oi^Oi^^ — 

If 11^=1 ''i / 1> then ai^^...^i^ = 0{ ,;^+^..+,J since n H h ife > n - sq 

"■1 

and wc assume that (2.2) holds for n — sq < k < n. By the above and that 
\ai\ < \/g(0),i > 2, then {ai^Oi^ ■ ■ ■ ai^^ai^+...+i^| = O(^) and therefore 

Mfe_i(t) =4ai+0 (^^j . 

Therefore, 

lim a'la^ = Mk-i. 

t—>oo 

By taking complex conjugate of the above, we obtain (2.2). □ 
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Furthermore, since (2.2) and that \ai\ < ^/g(0), i > 2, then lui^ai^ ■ ■ ■ \ 

0(\) and therefore 

a'la^-Mk-i = - iidnOia • • • aik(^h+-+ik = O ( \] ■ 

By taking complex conjugate of the above, we obtain (2.3). 

□ 

In the case that the lower Richardson moments vanish, better qualitative 
properties for the coefficients ai{t) and {ak{t)}k>2 are given as follows. 

Lemma 2.2. Let fn{^,t) = Y17=i ^ii'^)^^ ^ global strong degree n > 2 
polynomial solution to (1.1). Denote uq = min{A; > 1 | Mj. ^ 0}. Then 
(a) 

limaf'+Uno+i=Mno, (2.4) 

t—^0O 

lim a^afe = Mjfc_i, A; > no + 1, (2.5) 



t— »oo 



and 



(b) 



lim 

t—^0O 



lim a?°+^afe = 0, 2 < A; < no + 1. (2.6) 



a,{t) - V2t + Mo(0)] (2t)-o+3/2 = _^±1) |M„J2 . (2.7) 



Proof, (a) If no = 1, part (a) in Lemma 2.2 follows from Lemma 2.1. Now 
assume no > 2. We split the proof for no > 2 into two parts, 
claiml: 



-,no+l 
Hill C 

t—>^oo 

and 



lim 0"°+ Ono+l = , 
t—^oo 



lim o^afc = M/j_i, > no + 1. 

t— »oo 



Proof, (of claiml) As shown in (2.2). 



□ 



claim2: 



lim a?°+^afe = 0, 2 < A; < no. 
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Proof, (of claim2) In (2.1), the Richardson moments have the following rep- 
resentation: 

Mkit) = ^ iiai^ii^ ■ ■ ■ Oifc+iaii+.-.+ife+i 



E 



iiQi^ai^ ■ ■ ■ aii^^^ai^-\ — Hk+i- 



This means 



My 



E 



1'lO'hO'i2 ' ' ' 0'ik+i'^il-\ hife+i 



By taking complex conjugate, the above becomes 



1lO'ilO'i2 ' ' ' Oifc+ifliiH hifc+i 



■ (2.9) 



Hence 



O-k+l 



-.k+l 



Mk-Ol^ 



(2.10) 



By substituting no — 1 for k in (2.10) and applying the fact that M„q_i = 0, 
we can first show that 

lim a?°+^a„„ = 0. 

t^oo 

If no = 2, we are done. Assume no > 3 now and continue the proof by 
induction. 

Assume for some sq where < sq < no — 3, 



lim a^°+^a„o_s = 0, < s < 



t— >oo 



So. 



(2.11) 



We need to show that (2.11) holds for s = sq + 1. Denote Z = no — (so + 1). 
In (2.9), by substituting Z for /c + 1, we have that 



= — r 



E 



^1 Oil 0^12 ' ' ' ^ii^ii-i hi; 



(2.12) 
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However, if 11^=1 then ii + ■ ■ ■ + ii > Z + 1. Therefore, ai^+...+i; 

0( „o+i ) due to the assumption in (2.11), results (2.4) and (2.5). Hence 



iiajjOij • • • aj,aij+...+i, — a\ 



(2.13) 



Finally, by the representation (2.12) for ai and (2.13), we obtain the estimate 



ai 



a{-^0 



O 



and this implies 



lim a,a^°+^ = 0. 

t— »oo 



This says (2.11) also holds for s = sq + 1. 
Hence, claim2 is proven by induction. 



□ 



(b) 



By claiml and claim2, the proof for part (a) in Lemma 2.2 is done. 

ai{t) - VMo(O) + 2t 
al{t) - Mo(0) - 2t 



ai(t) + /Mo(0) + 2t 

-9{t) 



ai{t) + y/Mo{0) + 2t 
By (2.4), (2.5) and (2.6), we can get that 



lim9(i)(2t)"o+i = (no + 1) |M, 



no I 



Therefore, 



lim 



ai(t) - VMo(0) + 2i {2tro+3/2 ^ (no + 1) ^^^^p 



□ 
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2.2 Rescaling behaviors of global strong polynomial solu- 
tions 

Theorem 2.3. Let fn{^,i) be a global strong degree n>2 polynomial solu- 
tion to (1.1). Denote uq = min{A; > 1 | ^ 0}. We can see the following 
rescaling behaviors: 



(n) 



lim 

t—^0O 



lim 

t—^0O 



fniC,t) - VMo(O) + 2ie] (2*)"^ = M^,C'^\ (2.14) 

fn{V2t^,t)-(2t+^y-±mr,e 



{V2t) 



k=2 



Proof. Denote fn{^,t) = E7=i(^iit)C- 
(i) 

fn{C,t)-^/Mo{0)+2tC 

n 

= (ai{t) - VMo(0) + 2t) ^ + Y1 

i=2 

=1 + 11 

where I = (ai(t) - v^Mo(O) + 2t)^ and // = X;"=2 "-iC- 
By (2.4), (2.5) and (2.6) in Lemma 2.2, 

i+"o 



lim {2t) — iII)=Mn,e'^\ 
t— >oo 



By (2.7), we can obtain 



lim (2i)"o+i(/) =0. 

t—^0O 



By (2.16) and (2.17), we prove (i). 
(ii) 



fn{V2tc,t)=Y,aiit) (y^tjc- 



(2.15) 



(2.16) 



(2.17) 



i=l 
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^ ^ k=2 



k=2 



In order to prove (ii), it is enough to show the following claim. 
Claim: For A; > 2 



and 



lim 

t— »oo 



lim 



ai{t)V2i- (y2iy 



Mo(0) 



{2t) 



(2.18) 
(2.19) 



Proof, (of claim) The proof for (2.19) is easy. We will now focus on the 
proof of (2.18). Equation (2.8) states 



Mk-i — ^ iiOi^ai^ ■ ■ ■ aij.aji+...+jj, 
n.---.«fc;nLi»3^i 



We multiply the above identity by {y2t) and obtain 



ai 



Mk-i - J2 



ZlCljj^Qjj ■ ■ ■ fli^i-Oij-i 



By subtracting M^^i from the above identity, we have 

a^{V2if-Mk-i 



--Mk- 



hai^ai^ ■ ■ ■ a. 



0'i^O,i^-\ Hk 



ai 



^ ^ ilfljjGjj ■ ■ ■ 



(2.20) 
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In order to estimate the right-hand side of (2.20), we estimate 



2t 



a-i 



— 1 and 



in (a) and (b) as follows: 
(a) 

By the mean value theorem, there exists 9k{t) € [0, 1] such that 



2t 



ai 



-l = k 



2t 



k-l 



1 . 



For this term — 1, we have 



V2i ^ _ V2i -ai _ 1 2t - _ 1 g{t) - Mo(0) 

ai V2i + ai ai V2i + ai 



ai 



ai 



Therefore, 



1 = k 



\ ai 
This means 



2t 



ai 



- 1 



k-l 



(-Mo(0)+g(t)) 
al{^/2i + al) 



lim 

t—^oo 



2t 



ai 



-Mo(0) 



(k). 



(2.21) 



(b) 



E 



^lOil^jj ■ ■ ■ Oj^Clj-^-i 



= O 



1 



o 



1 



{2tf 



. (2.22) 



Now we come back to estimate the right-hand side of (2.20). By (2.21), (2.22) 
k 



\f2t 



- 1 



E 



l\(Xi^(Xi2 ■ ■ ■ OSi/jC'nH hife 



(2.23) 
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and 



E 



Zlfljjttjj ■ ■ ■ O^ij-fliiH Hk 



h,— ,H;nJ=i b'T^i 



2t 



(2.24) 



By (2.21), 



lim Mk-i 

t— »oo 



2t 



ai 



(2i) = -^^(A;)Mfe_i. (2.25) 



We apply (2.23), (2.24) and (2.25) to (2.20) and obtain that for k>2, 



lim 

t— >oo 



w(V2t) -Mfc_i 



Mo(0) 



ik)M, 



k-i- 



By taking complex conjugate of the above identity, we have for > 2 

k ^1 / Mo(0) 



lim 



ak(y2t^ -Mk-i (V2i) 



-kMi 



k-l- 



□ 
□ 



2.3 Geometric characterization of a rescaling behavior for 
moving domains 

Assume /n('^,i) is a global strong degree n > 2 polynomial solution to 
(1.1). We aim to give a geometric characterization for a rescaling behavior 
of Cl{t) = fn{D,t) as follows: 

Theorem 2.4. Assume fn{S,,t) is a global strong degree n > 2 polynomial 
solution to (1.1) with Richardson complex moments {Mj.}jt>o. Let Q,{t) = 
fn{D,t) and Q'{t) = {x|a/| ^{t) \ /irx € ^{t)} which has area tt. Denote 
no = min{A; > 1 | 7^ 0}. We show that 

1 I '^O 

lim max |b| - 11 (2t)^+— = |M„J 7^ 0, 

t^coz£dU'{t) 101/' 



1 I ^0 

lim max k(i, z) - 1| (2i)^+^ = (no - l)(no + 1) \M, 
t->ooze9n'(t) 



"■0 I 



where K(t,z) is the curvature for z G dQ'{t). This says that the decay rate 
^-(i+no/2) ^gg^ j„g^g ^g ^g^ ^j^^ ^^g j-^gQdiifig behaviors are precisely 

stated. 
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Proof. Denote fni^jt) = Y17=i^ii't)C- Let g{^,t) be the Riemann mapping 
where g : D ^ n'{t), g'{0,t) > and 5(0,*) = 0, then 



V2t + Mo(0)' 
(a)By (2.14), 



hm (2t)i+^ = M„„r+'- 

t— >oo 



Therefore, 

1 I ^0 

lim max - 1| (2t)^+^ = |Mno| / 0. 



(b)Let A{t) = ai{t) - A/2t + Mo(0), then 

g{^, t)=C+(f: , "'^'^ A + ( , ^^'^ i\ ■ (2.26) 
Here k — 1 can be expressed by g as follows: 



We set 

/n-l 

P=\ 



{k + l)ak+i k\ ^ I Mt) \ 
^2t + Mo(0) j l^72tTMo(0)/ 
(fc- l)(fe + l)afe+i A _ / ^(t) A 
^■^2 V2* + Mo(0) j \^/2t + m^) 



and obtain 5' (^, t) = 1+p and g'" (^, t)^ = p + q- Therefore, we express k—1 
by p, q as follows: 

1^1 V (1 + I5 I) / Iff I vi^ i(^ + i^ i) 1^1 / 

+ ^i?e - p4 + . . . ) + J_i^e {-qp + qp^ + ■■■). 
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By applying Lemma 2.2, we can obtain the precise decay rates of p and q to 
calculate the precise decay rate of — 1|. Here we treat two cases differently 
as follows: 

(Case i)Assume that Mi / 0. That is no = 1. The sharp decay rate of p is 

• If M2 / 0, then the sharp decay rate of q is f"^. Therefore the sharp decay 

rate of \k — 1| is 

• If M2 = and -M3 7^ 0, then the sharp decay rate of q is t~^l'^ . Therefore 
the sharp decay rate of |k — 1| is f^l'^ . 

• If M2 = 0, M3 = and M4 / 0, the sharp decay rate of q is . Therefore, 

the sharp decay rate of |k — 1| is 

• Finally, if M2 = -M3 = M4 = 0, then the sharp decay rate of q is t~^. In 
this case, the sharp decay rate of |k — 1| is t~^. 

In (Case i), it is clear that 

lim max k(t, z) - 1| (2t)^+ 2 = = (no - l)(no + 1) |M„„ I . 
t^oo 2gaQ'(t) 

(Case h) Assume Mi = 0. That is no = min{A: > 1 | M^ 7^ 0} > 2. The 
sharp decay rates of q and p are both ^^(^+"0/2) jj^ ^j^jg c^se, the sharp 
decay rate of |k — 1| is 

k-ll = 



2 



Moreover, in (Case ii). 



1 J- ^0 

lim max kft, z) — 1| (2t)^ 2 
= lim ma^|i?e(g(C,t))|(2i)^+^ 



■- lim 
t— +00 



(no - l)(no + l)ano+i 



V2t + Mo(0) 
(no-l)(no + l)|M„J. 



(2i) 



1+^ 



□ 



3 Rescaling behavior in Stokes flow 

Assume Y17=i ^j(^)?* ^ global polynomial solution to zero surface tension 
Stoke flow driven by injection with strength 27r. The invariant quantities 
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Cfc = Yl]=i j^j^j+k for 1 < A; < n — 1 as shown in Cummings, Howison and 
King [1]. Here Co{t) = J2]=ijaj{t)a~(t) = Co(0) + 2t. Since \aj{t)\ < jai{t), 
then 

n n 

Co(0) + 2t = J2jaj{t)^ < J] 

i=i i=i 

which impUes that af(i) > .3 (2t + Co(0)). 

3.1 Decay of coefRcients in terms of invariant quantities 

Lemma 3.1. If fn{£,-,t) = Z^ILi ^iC* is a global degree n> 2 polynomial 
solution to (1.2), then ak{t)ai{t), k > 2 is bounded. 

Proof. It is sufficient to prove that On-jOi is bounded for < j < n — 2. 

Since C^-i = aio^. Therefore the above result holds for j = 0. Now 
assume that a„_jai is bounded for < j < sq- 

Since 

Cn-(so+2) = «.l«'l+(n-so-2) + ^ 30'j<^j+{n-so-2)^ 

3=2 

therefore, 

so+2 

ai«n-(so+l) = C'n-(so+2) " X] jaj«j+(n-so-2) ■ (3-1) 

J=2 



Here, for 2 < j < sq + 2, 

/ ao\ , s I I 

< j |aiaj+(„_so_2)| 



ri%+(n-so-2)| 



which is bounded since we assume an-jCLi is bounded for < j < sq. There- 
fore, Yli^j^=2 j^j(^j+{n-so-2) is bounded. We obtain that aian-(so+i) (3-1) 
is bounded. □ 

By (3.1) and similar argument as that of Hele-Shaw flows, we obtain the 
following results. 

Lemma 3.2. If fn{C,t) = Y17=i^i^^ ^ global degree n> 2 polynomial 
solution to (1.2), then for 2 < k < n, 

ri-(fc-l) 

^Ihn [aj^ai - Ck-i] {V^if = ^ jCj-iCj+fe-2 
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and 




k-l 



2 



k=2 

Theorem 3.3. If fn{^,t) = J2^=i^i^^ ^ global degree n>2 polynomial 
solution to (1-2), then 



3.2 Geometric characterization of a rescaling behavior for 
moving domains 

Given fniS,,t) which is a global degree n > 2 polynomial solution to (1.2). 
We aim to give a geometric characterization for a rescaling behavior of 
^{t) = fn{D,t) as follows: 

Theorem 3.4. Assume fn{^,t) is a global degree n> 2 polynomial solution 
to (1.2) with invariant quantities {Cfc}fc>i. Let fl{t) = fn{D,t) and Q'{t) = 
{x\^/\Tl(t)~\~/^^x G ^{t)} which has area tt. Then 



n 



^2t+c^m V2t = Y,Ck-i(!'- 



k=2 



n 



lim max 1 1 z 
t-^co zedfl' (t) 



1| (2t) 



max 



k=2 



lim max 1^(^,2;) 




where K{t, z) is the curvature for z G dW 



it). 



Proof. The proof is similar to that of Theorem 2.4. 



□ 



16 



Acknowledgements 



This paper is an extraction from my Ph. D thesis. I am indebted to my 
advisor, Govind Menon, for many things, including his constant guidance 
and important opinions. This material is based upon work supported by 
the National Science Foundation under grant nos. DMS 06-05006 and DMS 
07-48482. 

References 

[1] L. J. CUMMINGS, S. D. HowiSON, AND J. R. KiNG, Two-dimensional 
Stokes and Hele- Shaw flows with free surfaces, European J. Appl. Math., 
10 (1999), pp. 635-680. Hele-Shaw flows and related problems (Oxford, 
1998). 

[2] B. GUSTAFSSON, On a differential equation arising in a Hele-Shaw flow 
moving boundary problem, Ark. Mat., 22 (1984), pp. 251-268. 

[3] B. GusTAFSSON, D. Prokhorov, and a. Vasil'ev, Infinite lifetime 
for the starlike dynamics in Hele-Shaw cells, Proc. Amer. Math. Soc, 
132 (2004), pp. 2661-2669 (electronic). 

[4] C. HUNTINGFORD, An exact solution to the one-phase zero-surface- 
tension Hele-Shaw free-boundary problem, Comput. Math. Appl., 29 
(1995), pp. 45-50. 

[5] O. S. Kuznetsova, On polynomial solutions of the Hele-Shaw problem, 
Sibirsk. Mat. Zh., 42 (2001), pp. 1084-1093, iii. 

[6] S. Richardson, Hele-Shaw flows with a free boundary produced by the 
injection of fluid into a narrow channel, J. Fluid Mech., 56 (1972), 
pp. 609-618. 



17 



